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Abstract: We present a mesoscopic model for thermoelectric phenomena in terms 
of an interacting particle system, a lattice electron gas dynamics that is a suitable 
extension of the standard simple exclusion process. We concentrate on electronic heat 
and charge transport in different but connected metallic substances. The electrons 
hop between energy-cells located alongside the spatial extension of the metal wire. 
When changing energy level, the system exchanges energy with the environment. 
At equilibrium the distribution satisfies the Fermi-Dirac occupation-law. Installing 
different temperatures at two connections induces an electromotive force (Seebeck 
effect) and upon forcing an electric current, an additional heat flow is produced at 
the junctions (Peltier heat). We derive the linear response behavior relating the 
Seebeck and Peltier coefficients as an application of Onsager reciprocity. We also 
indicate the higher order corrections. The entropy production is characterized as the 
anti-symmetric part under time-reversal of the space-time Lagrangian. 

KEY WORDS: interacting particle system, Seebeck and Peltier effect, fluctuations 
and response functions. 



1. Introduction 

The goal and motivation of the paper comes form exploring nonequilib- 
rium behavior via the theory of spatially extended stochastic dynamics, 
a standard reference is ^J. We present an interacting particle system 
for the standard thermoelectric phenomena such as the Seebeck and the 
Peltier effect. Second, we show how the steady state fluctuation theorem 
rigorously reproduces the phenomenological linear response theory. Fi- 
nally, the model also provides a testing ground for deriving higher-order 
response functions and for exploring the nonequilibrium theory beyond 
linearity. 



Thermoelectric systems provide a phenomenology for the foundations of 
nonequilibrium thermodynamics. Standard examples of thermoelectric 
phenomena are the Peltier effect (1834) — electric current through differ- 
ent metallic conductors will in general cause production or absorption of 
heat at the junctions, and the Seebeck effect (1826) — an e.m.f. (electro- 
motoric or thermoelectric force) will appear in the circuit if the junctions 
are maintained at a different temperature (thermocouple). They are 
discussed in the standard textbooks such as [2] and except for a short re- 
minder in Section we assume the reader familiar with the basic facts. 
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As written in the first line of Onsager's 1931 article on reciprocal rela- 
tions, "When two or more irreversible transport processes (heat conduc- 
tion, electrical conduction and diffusion) take place simultaneously in 
a thermodynamic system the processes may interfere with each other." 
Reciprocity relations are the most well-know realizations of that and 
Thomson was probably the first (in 1854) to propose them. Onsager 
took up the derivation of the reciprocity relations now named after him 
[H]. In the second half of the previous century the linear response the- 
ory added the response relations such as these of Green and Kubo from 
which the Onsager relations are explicitly visible. 

The more recent history has seen a revival in attacking questions on the 
construction of nonequilibrium statistical mechanics. One proposed gen- 
eralization of (linear) fluctuation-dissipation relations has been the so 
called fluctuation theorem, both in its transient and in its steady state 
version. For steady states, it gives a symmetry in the distribution of the 
variable entropy production. We derive it here in Section 15.21 Section |H1 
explains how it leads to the phenomenological relations between Seebeck 
and Peltier coefficients. Not emphasized in the literature so far is that 
this fluctuation symmetry is useless for higher order response functions; 
it is unable to produce a useful relation for the second order derivatives of 
the current around equilibrium. We observe in Section [7| how to remedy 
the situation for our particular model. 

The physical origin of the model is discussed in Section 13.11 From the 
point of view of the standard theory of interacting particle systems, it 
is an adaptation of exclusion processes, see [4^ and i5j for mathemati- 
cal background. In that way it is easy to simulate and it thus offers an 
environment for studying analytically less accessible aspects of thermo- 
electric phenomena (such as the influence of disorder) . From the point of 
view of solid state physics, it can be seen as the result of time- dependent 
perturbation theory to transport problems in metals. We only keep the 
electron hopping with transition probabilities that satisfy a local detailed 
balance condition. 

The details of the model are explained in Section 01 while, as a warm-up. 
Section |21 contains its most simplified form. 

2. Toy model for heat conduction 

Consider a lattice rectangle of height 2 and width M, see Fig.l. Each 
of the 2M sites can be occupied by at most one particle. Think of 
the rectangle as a 'metal wire' of size M in which each couple of sites 
(x, 2M -|- 1 — x) forms one cell with two energy levels that can each be 
occupied by one electron. The energy difference between the two levels 
is taken equal to a constant e > 0. 

We consider a Markov process on the state space K = {0, Ij^A/ which 
particles can hop on nearest neighbor sites. The hopping is symmet- 
ric (subject only to the simple exclusion rule of one particle per site) 
except for the bonds (2M, 1) and (M, M -|- 1). More specifically, an 
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elementary transition consists of randomly picking a site x and chang- 
ing the configuration r] into the new configuration T^r] which remains 
equal to r] on all sites y except for y = x,Txrj{x) — r]{x + 1) and for 
y — x-\-l, Txr]{x + 1) — r]{x). The transition rate for rj — > T^r] is equal to 

ex.p[/3Ae{7]{x + l)-r]{x))/2], if x = 2M 

ex.p[(3B e {r){x) - r){x + 1)) /2], if x = M 

1, otherwise (2.1) 

The number of particles is conserved in the evolution and the Pa and 
Pb are the only (other) parameters in the model. In words, the particles 
jump to nearest neighbor cells but remain on the same energy level (no 
effective scattering in the bulk). In the left and right end cells, particles 
undergo thermal transitions; they can change energy level at inverse tem- 
peratures Pa and Pb respectively. The system is thus at its both ends 
in contact with a heat reservoir and one expects a heat current flowing 
whenever they have unequal temperatures. 

[FIGURE 1] 

The model satisfies the condition of detailed balance when left and right 
Pa = Pb = P- For each function F of the total particle number, the 
probability 

^ 2M M 

Puiv) = -^^Q^vix)) exp[-pe^ri{x)], rj e K 

x=l x=l 

is stationary and reversible for the dynamics. These equilibrium weights 
are completely determined by the number of particles X^fli v{^) the 
upper energy level. In order to make our calculations possible we make 
use of the equivalence of ensembles to express the single site marginal in 
terms of chemical potentials. In equilibrium the grand canonical setup is 
made by attaching particle reservoirs with chemical potential fJ-A — fJ-B — 
fi to the left and right side of the system. For this setup we now make 
the choice F{n) = exp[p /in]; is a homogeneous product measure with 
single site marginal 

n,ix) ^ Prob,, [rjix) = 1] = ^^^^J^^^ ^ ^ (2.2) 

where e{x) = or e depending on whether the x is 'down' M + 1 < x < 
2M or 'up' 1 < X < M in the energy cell. 

If the reservoir temperatures are unequal, when Pa 7^ Pb, there is no ex- 
plicit expression for the stationary measures but we do know the density 
profile n{x) = PTohp[ri{x) = 1] in a stationary measure p. In the bulk 
and because of the symmetric exclusion rule, there is a linear density 
profile in both the upper and the lower energy layer: 

^(^) = ^u + ^f^K-n^] l<x<M 
2M ~ X 
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but the density in each cell (containing the couple (x, 2M + 1 — x), 1 < 
X < M) is constant for a canonical system due to particle conservation: 
n(x) + n{2M + 1 — x) = + = + n^- In the steady state, there 
is no particle current between cells but there may be a heat current. 

On the upper energy layer particles transport an energy equal to e; giving 
rise to a steady energy current, taken negative in the x— direction since 
we will heat up side B, equal to 

R A 

W = e^— f (2.3) 

Remember that and are steady densities (left and right on the 
upper energy layer depending on M,Pa and Pb- On the other hand, 
as M I +00 (very long wire) these densities tend to their respective 
equilibrium values as for example from ()2.2|) . 

I I (2 4) 

exp(/?^[e-/iA]) + l " exp{pB[e - IXb]) + I 

That is an application of convergence results for a sequence of Markov 
processes, see e.g. Chapter 4.9 in [H]. It should be emphasized that fj,A 
and /iB are a function of Pa and Pb respectively, and of course as well 
of the particle density. An expression for fx a and fiB may be derived by 
solving the following identity: 



exp(/5ij[e - ij,b]) + 1 exp(/5B[-yUB]) + 1 
Naturally the identity for side A is similar. 

As a result of ()2.4|) . the heat current ()2.3|) satisfies the stationary Fourier 
law with 

< K = limM(jQ) < +00 (2.5) 

Writing out the right-hand side of ()2.3|1 with the substitution ()2.4j) for 
Pa = 1/7" and Pb = 1/(T + e), we get a heat conductivity equal to 



with H{e) = 1/4 cosh^(/?(/i — e)). The last term is the contribution due 
to the canonical ensemble whereas the first two terms would also occur 
in a derivation for a grand canonical setup. 

We conclude that this toy-model already shows electronic heat conduc- 
tivity in a clear and simple way. 



The next Section extends the model by allowing multiple upper and lower 
energy levels and by interpreting the upper and lower levels as character- 
izing two different but connected metal wires, see Fig.2. Again at the two 
junctions there may be different temperatures and we will also consider 
the case where electric power is applied. 
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3. The general model 

3.1. Physical input. Metals are excellent conductors of both heat and 
electricity. Already the empirical law of Wiedemann and Franz (1835) 
relating thermal and electrical conductivities, points to the electrons as 
heat carriers. While say for the specific heat of a metal at room temper- 
ature there is no observable electronic contribution, for nonequilibrium 
purposes mostly electrons are responsible for heat and charge transport. 
As the electron mass is small, even at room temperatures and for typi- 
cal metallic electron densities, the Maxwell-Boltzmann and Fermi-Dirac 
distributions are very different. Therefore, the Pauli exclusion principle 
must play a role and the quantization into energy levels, even if only 
near the Fermi energy, must be part of any physical model of transport 
in metals. 

On the other hand a fully quantum mechanical and microscopic treat- 
ment of transport problems is typically out of reach. We have here a 
system with extremely many degrees of freedom and one usually finds 
refuge in perturbation theory. The electrons are moving almost indepen- 
dently but subject to (weak) interactions with phonons, impurities and 
imperfections. When the average time between collisions multiplied with 
the Fermi energy is larger than Planck's constant (as is the case in met- 
als), Fermi's Golden Rule can be applied to the transition probabilities 
between electronic states and in ways similar to that for the Boltzmann 
equation, one obtains a Master Equation. 

The resulting model that we take up in the present paper is mesoscopic 
(only electrons moving on a lattice) and while it is thus physically moti- 
vated from time- dependent perturbation theory for quantum mechanical 
processes, it fits perfectly well in the classical theory of interacting par- 
ticle systems. 

One imagines two wires connected at their ends. The wires (long thin 
bars) consist of a large array of cells, discrete lattice points, at which 
are located a finite number of energy levels. Free electrons are hopping 
between nearest neighbor cells, subject to an exclusion rule and with 
transition rates that are specified via a local detailed balance condition. 
We start by introducing the general set-up. The simplest version was 
already in the previous section and other simplifications follow later, see 
Section Eini 



3.2. State space. Consider a ring of N points or cells. The situation 
one should keep in mind is that of connecting two different wires, the 
first of length M, the second of length — M. To each cell in the 
first wire, 1 < k < M, is assigned a set Si of > 1 energy levels 
, and to each cell in the second wire, M + 1 < A; < A^, is assigned 
a set S2 of ^2 > 1 energy levels. The set of energy levels at cell k 
is in general denoted by £k but it should be understood that it equals 
Si if k e {1,...,M} and that St = S2 for k e {M + 1, . . . , N}. A 



6 



general element of is written as e, e', These are the electronic 

states (we ignore spin). There can be at most one particle per energy 
level and we write x = {k,e),r]{x) = r]{k,e) G {0,1} for e G Sk- The 
total configuration space is thus K — x^^^jO, 1}'^^ >^k=M+i {0> ^}^^ 
its elements are denoted by r],^, 

[FIGURE 2] 

3.3. Transitions. In each of the two wires particles can hop within the 
same energy level to nearest neighbor cells. That is subject to the ex- 
clusion rule and possibly biased in one direction via some external field. 
We write > for the electric field over the bond {k, k + 1). One can 
think of a gradient in chemical potential which is locally installed but 
the details of the "battery" will be of no concern here. One imagines 
that the two wires are in a heat bath at inverse temperature (3 (Boltz- 
man's constant is set equal to one) except possibly for the junction 
(M, M + 1) which is kept at a different temperature. 
Fixing A; = 1, 2, ... , M-1, M+1, . . . , for the "bulk" bonds {k, k+1), 
the exchange operator associated to energy level e G £fe is 

{V,^rj){k + l,e)^rjik,e), 

{V^r]){x) — r]{x) otherwise 

and such an exchange of occupations 

T] — * V^T] takes place at rate exp[/3 Ek{ri{k, e) — ri{k + 1, e))/2] 

At the junctions, k = M or k = N, particles can also hop but now to an 
arbitrary energy level of the other wire. These transitions are described 
via the exchange operator T^'^ for e & Ek,e' & £k+i- 

(Tr'^)(A:,e)=^(A; + l,e'), 
{T/v){k + l.e')=r^{k,e), 

(T^'^ r]){x) — r]{x) otherwise 

and 

77 ^ Tl'^'i] at rate exp[/?^'(^fc + e - e'){r]{k, e) - r]{k + 1, e'))/2] 

describes a local detailed balance condition. We set (5^ = (3 and at the 
other junction — P + A. A 7^ means that one junction is main- 
tained at a different temperature. 

Obviously the Vk are just a special case of the T^. An immediate gener- 
alization would be to allow transitions with T^'^ everywhere, also in the 
bulk, or to allow on-site thermal exchanges. That would break bulk en- 
ergy conservation but it would allow to describe thermal exchanges with 
the environment at all cells, e.g. from scattering. In the same sense, it is 
no problem to add more jimctions (connecting more than two wires). As 
we have it now, the bath's role in the bulk is purely as a work reservoir. 
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One could also imagine adding particle reservoirs. By that we mean se- 
lecting one or more cells k on which at energy level e G we have the 
transitions 

{Alri)ik,e) = l-r]{k,e), 

{Alr]){x) = rj{x) otherwise 

These break particle conservation: particles flow in or out with rates 

r] ^ Air] at rate a^. exp[/3(e — /ifc)(?7(/c, e) — |)] (3.1) 

where /ifc is the chemical potential. We will not deal with these gener- 
alizations more explicitly. It is possible to include them but soon the 
notation starts to get heavy. Let it be clear however that in no way do 
we need to introduce a temperature or chemical potential associated to 
the internal local conditions of the cells in the wire. Such close to equi- 
librium conditions are mostly present in phenomenological treatments of 
thermo-electric effect. 



3.4. Dynamics. We consider a continuous time Markov process on K 
consisting of the two types of elementary transitions T^'^ and intro- 
duced above. For fixed sets of energy levels Sk and spatial extensions M 
and N, the remaining parameters are (electric field at bond {k, k + 1)), 
P and A (inverse temperatures). 

For Ek > work is being done on the system and the particles jump 
preferentially forward, i.e., in the direction k —>■ k + 1. Energy is not 
(necessarily) preserved in the transition T^'*^ but the total number of 
particles is conserved (at least when the = in ()3.1|) . which we will 
assume hereafter). 

The backward generator L of the Markov process is a sum 

L= J2 Lu + Lm + Ln (3.2) 

with 

Lkfiv) = E ^Mi^EMK e) - vik + 1, e))/2] [/(V^^r?) - /(r/)] 

LNfiv) = E exp[/?(E^+e-e')(r7(iV,e)-r/(l,e'))/2][/(T^'^'r/)-/(r/)] 
and 

LmKv) = E exp[{(3+A){EM+e-e'){r]{M,e)-r]{M+l,e'))/2][f{nfr])-f{r])] 

ee£M,e'e£M+i 

Denoting the general rate as r{ri, rj') (of course depending on all possible 
parameters), we have 

Lf{r]) = J2^{v,v')[f{v')-f{v)] (3.3) 

V' 
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summarizing ()3.2|) ; the Master Equation for the probabihties at time t is 
as usual 

^Probt(?7) = X] ['^('^''^) Probt(?7') - r{ri,r]') Probt(r7)] 

ri' 

For every initial configuration rjQ the above defines the process rit,t > 0. 
We are interested in a stationary regime for fixed (but imagined very 
large) spatial extensions M and N. Let p denote a stationary measure 
for the above Markov process. It could depend for example on the particle 
density. 

We write for the stationary process over the time-interval [0, r]; uj = 
{.Vt)t=o is a path as sampled from P^. We will assume that there is a 
finite positive constant C with 

for each pair rj, rj' that can be connected via the transitions outlines above. 
In general we know very little about this stationary process; under equi- 
librium conditions we can make it more explicit. 



3.5. Equilibrium. The system is in equilibrium when A = and all 
Ek = 0. Consider the probability measure p/j,^ defined as 

1 ^ 
PpAv) = -e^'^S^''(^)exp[-/3^^er/(A;,e)] (3.5) 

fc=l eGSf^ 

That is a product measure but it is not (in general) homogeneous. It is 
a stationary reversible measure for the equilibrium dynamics. Its single 
cell marginal is given by the Fermi-Dirac distribution, for e & Sk, 

Pp,Mk, e) = 1) = ^^(,_^) ^ ^ (3.6) 

That corresponds to the grand-canonical set-up. Except when we add 
sources and sinks for the particles (a^ 7^ 0), the total number of particles 
is conserved in the dynamics as defined above. 

Being in steady equilibrium implies ignoring the difference between ther- 
modynamic past and future and on average there is no net transport nor 
dissipation. 



3.6. Simplification. The simplest version is obtained by taking £^ = 
{ci} and = {62} singletons. There is then one energy level per cell, 
rii = n2 = 1, and we may just as well put 62 — 61 = e. We also put 
ak = = Ek for all k (canonical set-up without external field). In the 
bulk of each of the wires we now only have symmetric hopping and at the 
two ends of the wires (we are now talking about the bonds (M, M + 1) 
and {N, 1)), there are thermal transitions: 

t]^TmV at rate exp[{(3 + A) e {7]{M) - r]{M + l))/2] 
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through which the occupations are exchanged over the bond (M, M + 1) 
and 



through which occupations are exchanged over the bond (A^, 1). We see 
that it formally corresponds to the toy model for heat transport of Section 
121 with 2M = N and = /5, /3b = /5 + A at the junctions. In that way, 
the model system of Section |21 has two possible physical interpretations, 
one as consisting of M cells each with two energy levels, the other as a 
couple of two wires each of length M and each with one energy level per 
cell. In that last interpretation, the heat current in Fourier's law (defined 
in ()2.3p ) is in fact a particle current through the upper wire, caused by 
the different temperatures Pa 7^ Pb at the junctions and so provides the 
simplest illustration of a Seebeck effect. To see the Peltier effect, we take 
A = (equal bath temperatures) but some > 0. A particle current 
arises through the ring and energy of magnitude e gets dissipated at the 
junctions. 



4.1. Local particle and heat currents. The basic physical quantities 
in the phenomenon are energy and particle number. The electric or par- 
ticle current consists of particles hopping to nearest neighbor sites. We 
fix a path u = {r]Q, . . . ,77^) in which each change (at a random time) 
corresponds to one of the transitions described under Section 13.31 
The net particle current (integrated over a time-interval [0,r]) from en- 
ergy level e at cell k to energy level e' at cell k + 1 equals 



for e G Sk,e' G Sk+i- The sum in the right-hand side of (j4.H) is over 
all jump times t; the moments at which the trajectory uj changes states. 
Here and in what follows, all times in the sums are t = t~ , right before 
the jump. So ()4.1|) is a sum of +l's and — I's according to whether the 
particle jumps from {k, e) to {k+ 1, e') or oppositely; the jump is counted 
at the time of the transition. When the bond {k, A; + 1) is in the bulk of 
the wires {k 7^ M,N), then only transitions between levels of the same 
energy are possible and the net number of particles jumping forward 
equals 



T] — » Tj^T] at rate exp[— /9e {ri{N) — ri{l))/2] 



4. Currents and dissipation 



J\e,e'){u;)^J2^l^t{k,e)~r]t{k + l,e' 



')] 



(4.1) 



J'{e, e)(u;) = r^'iu) ^ ^[r/,(fc, e) - r^,{k + 1, e)] 



t 



We call 





(4.2) 



the variable electric current over the bond (fc, k + 1). 

Secondly, energy is transported. Over the bond {k, k + 1) in the bulk the 
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time-integrated energy current equals 

At each junction k = M, N, there is a local heat exchange with the 
environment, commonly called the Peltier heat (counted positive for a 
heat transfer into the environment), here variable and equal to 

4^ J2 {e-e')j\e,e') 



4.2. Conservation laws. In a configuration r], 

are the energy, respectively the particle number at cell k. These change. 
We fix again a trajectory u = (r/o, • • • , Vr)- For k — 1, k ^ M, N 

UkiVr) - Ukivo) = Jh~\^) - 4(^) (4.3) 
At the junctions k = M, N, 

Uk{vr)+Uk+i{vr)-Uk{m)-Uk+i{vo) = 4''(^)-4^'(^)-4M (4-4) 

Similar relations hold for particle conservation. In particular, 

NkiVr) - Nk{r]o) = Jlr' - Je (4.5) 
As a consequence, the steady state average defines the particle current 

Je ^ (4) 

independent of k. 

In the steady state (j4.4j) can be rewritten as 

(4"') = (4) + (4""') 

with, from (j4.3p . {J^) = Ji constant for /c = 1, . . . , M - 1 and ( J^) = J2 
constant for fc = M + l,...,A^ — 1. Therefore, denoting the Peltier heat 
at junction (M, M + 1) with Jb = {Jp) and at junction {N, 1) with 
J4 = {Jp), we have 

Ji = Jb + J2, Ja + Jb = (4.6) 



4.3. Heat dissipation. The first law of thermodynamics 

dU = -dQ + dW (4.7) 

with dQ the dissipated heat and dW the work done on the system can 
be applied to our model by identifying the proper currents. The total 
energy of the system 

N 
k=l 

is not globally conserved, see ()4.4|) . The Joule heat is the (electrical) 
work done 
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over the bond {k, k + 1). As a result, for a path u = (?7o, . . . ,r7r) the 
variable dissipated heat over a time-interval [0, r] is 

N 

k=l k=M,N 

Since each bond {k, /c + 1) is assumed attached to a large heat bath, we 
get that as a function of the system's history the total change of entropy 
in the environment equals 

N 

S{u) = J2 (^Ek4{u) + 13 {u) + (/3 + A) (4.8) 
fc=i 



5. Entropy production 



5.1. Driving forces and mean entropy production rate. There are 
in fact two driving sources of nonequilibrium: one is the temperature 
gradient, the others the electric field. The temperature gradient is in- 
stalled by putting different temperatures at the junctions. The electric 
field models the presence of an external electric field or more simply, of 
a battery and the installation of an electric potential. 
In steady state average, the steady entropy current equals the mean en- 
tropy production 

{S)=[3Y, Ek{J%) + (/? + A) Jb + (3 J A (5.1) 

A; 

as obtained from ()4.8|) . Writing 

and using the conservation law ()4.(jj) . we have 

1^{S) = P[QJe-VTJs\ (5.3) 

for the entropy source strength per unit length of the wire. The driving 
forces are clearly visible: fi is the gradient of electric potential and VT 
is the temperature gradient. J5 = (/3 + A) is the entropy flux into the 
reservoir controlling the junction at inverse temperature /3 + A. 

The mean entropy production rate a over a time-interval [0, r] is defined 
as 

I f . . 1 

= a = - 



S{uj)dP;{uj) = -{S) (5.4) 

T J ' T 

and can of course be rewritten from ()5.3p . We have the fundamental 

result that it is positive: When not in equilibrium cr > 0. That can be 
seen as follows. We compute 

fl-log^ (5.5) 
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where G is the dynamical time-reversah Bo; = {'r]r-t)t=o- Then, by 
construction, 

J (iP;(cu)e-^(^) = 1 

and hence, 

J dPl{uj)R{uj) > (5.6) 

and is non-zero whenever R{uj) is not constant (with PJ-probabihty one). 
The hypothesis "not in equihbrium" refers exactly to that breaking of 
time-reversal invariance. The computation of R is an application of the 
Girsanov formula for Markov jump processes, see e.g. the Appendix in 
|3]. It is easy to check from Section 3.4 that 

R{uj) = S{uj) - \ogp{r]^) + \ogp{r]o) (5.7) 

and hence the steady state expectations of R and of S coincide. As a 
result, (|5.6|) implies the positivity of (|5.4|) . 

5.2. Fluctuations. Consider again ()4.8|) : 

S = S{E, A) = So + J2 ^'^'^E + ^Jp, So = (3{J^ + J^') (5.8) 

k 

In equilibrium = A = and always {Sq) = 0. The dissipation ()5.8|) 
is a fluctuating quantity and we can consider its generating function 

Clearly 5*0 = — S* so that by combining with (j5.5p and (j5.7|) . we get 

From assumption (j3.4p and for real z that implies the fluctuation sym- 
metry 

p{z)=p{l — z), p{z) = — lim —\ogFr{z) (5.9) 

We can however look in somewhat greater detail to the fluctuations as 
they arise from the local currents. 

Given numbers 7'^(e, e'), we write 

N 

k=l e€£k,e'&£k+l 

For example, with 'y''{e,e') = ip^{e,e') = Ek + {e — e')l3'^, ■ J = S 
the entropy current. If we take 'y^{e,e') = {pQ{e,e') = (e — e')/3, then 
ipo ■ J = So the equilibrium value defined in ()5.8|) . 

Consider now the following generating function for the fluctuations of 
the currents 




(5.10) 
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The second argument of p (the dependence on Lp) comes of course from 
the dependence of the steady state law = dP^ on the driving parameters. 
The hmit exists for the relevant choices of (7'^(e — e')) via an application 
of the Perron- Frobenius theorem, see e.g. jH]. Moreover, 

p(7, ^) = p(^ - (^0 - 7, (5.11) 

To see it, we start again from the definition ()5.5|) for R and observe that 
for all functions /, 

j dP;{u)f{Qu) = J rfP;(c.)e-^(-) fiuj) (5.12) 

We substitute / = exp 7- J and use first that J''(e, e')(9a;) = — J^(e, e')(to'), 
the antisymmetry under time-reversal of the variable currents. On the 
other hand, R = S — log p(?7r) + log ^(^70) with 

S = {(^- ipo) ■ J + So 

given in ()4.8|) or in ()5.8|) where, as follows from ()4.3|) - ()4.4|) . the equilib- 
rium value 

Soiuj) = 00 ■ J{oo) = (3 [U{r]^)) - U{r]o)] 
equals the total change of energy. 
For that choice of / 

C j dP^e-^^-^^'^-^f] < j dP;e-^f (5.13) 

where we used again ()3.4|) . We conclude that after taking the logarithm 
and dividing by r t +cx), R can be substituted with {0— (po) ■ J in (j5.12p 
and the identity ()5.1ip follows immediately. 



6. Linear response regime 



6.1. Green-Kubo relations. The steady state average will from now 
on be decorated by the field (p to indicate its dependence on the driving: 
(■) = (■)>^- For the equilibrium values (po we write = (■)o. A 

general couple of energy values e G £k,e' G £k+i is denoted by a = 
(e, e'), J^{e, e') = J^, (p''{e, e') = (p''{a). In particular, for all a and k, 

{j')o = 

is the zeroth order contribution to the steady currents. The Green-Kubo 
relations want to give details about the linear order in around 0o'- 

s,e 

up to higher order terms, with 

g A: 

La,k;5/ = {J a) 0=00 (6-1) 

linear transport coefficients. 
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The following Green-Kubo relation follows directly from ()5.11|) : 

L^,kM = h lim -(4'4)o (6.2) 

rt+oo T 

where one should remember that the currents are extensive in r. 

Deriving ()fi.2j) proceeds by differentiating the left- and right-hand sides 
of (jnHH). With the notation 

C^<?(7,^)^ A^(7,^) (6.3) 



^'^a^l-^^^?(7,^l (6-4) 

and under the condition that the function p in ()5.10|) is smooth, we have 
in fact 

Z}^C^p(0,0) = -iC,^CMO,0) (6.5) 

which implies (16. 2j) . The reasoning above is a detailed reproduction 
of what was observed before, e.g. in [3 IHl Ej, to derive fluctuation- 
dissipation relations from the fluctuation symmetry (jS.llj) . One observes 
that the Green-Kubo relations ()6.2p immediately imply the Onsager reci- 
procity 

La,k;5,e = Ls,e;a,k (6.6) 

6.2. Reproducing the phenomenological behavior. The above Green- 
Kubo relations ()6.2|) are the most detailed forms of linear response rela- 
tions. We connect them now with the standard experimental set-up for 
the situation of the different thermo-electric effects. 



The basic equation is (j5.3p . Though we did not start from a balance 
equation for the "close to equilibrium entropy of the system," (j5.3p re- 
produces exactly equation (61) in Chapter XIII of |2]. The next thing is 
to linearize the currents Js and Je around VT =6 = 0. For that we have 
the Green-Kubo relations. We can indeed either work from the previous 
section or, more directly, use ()5.12|) which says that 

(J) = i(J(l-e-^)) (6.7) 

for functions J that are antisymmetric under time-reversal, J0 = — J, 
like currents. 

We specify to the case where just one of the bonds, pick {k, k + 1), carries 
an electric field Ej. = EN (proportional to the length of the wire), all 
others being exactly zero. That corresponds to the situation where locally 
in the wire some external electric power is applied to particles that do 
not interact besides exclusion. This may seem unrealistic but it keeps 
things manageable and we expect the results to be equivalent to a more 
realistic description. We then have, up to linear order in E and A, 

. u,„i[^?(/!±^(,«4>„ . ^^m] (6.8) 

T T T T z z 
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Rewriting that for js = lini^ Js/Nt (per unit length and per unit time) 

js = -Lu VT + Li2 E (6.9) 

with 

Lu^ f hmi((40^)o 
/ T r 

L^2^ ^^^±^hmi(J*V|)o (6.10) 

In the same way, 

hm i = hm 1 [^(( J^)^)o + ^{J%J^)o\ (6.11) 
T r T r Z Z 

so that for je = hm,- Je/Nt, 

je = -L2iVT + L22E (6.12) 

with 

^22= flim-((J|)^ 
2 T r 

L2i^ ^^^^lim-(J*^j|)o (6.13) 
2 T r 

The Onsager relation L12 = L21 is apparent and hence, if we write the 
standard form 

Js = -\PVT + 71PJE 

E = -vVT + RjE (6.14) 

we get A/3 = Lu — L21L12/ L22 for heat conductivity X, R = I/L22 is a 
resistivity, and the coefficients TrjS = L12/L22 and v = —L21/L22 verify 
the second Thomson relation 

V = —7T (3 

The coefficient v is called the differential thermoelectric power and 7r/3 
is the entropy current per electric current when temperatures are equal, 
see 12]. 

6.2.1. Seebeck effect. That is the effect of so called thermoelectric power 
that can be measured via a thermocouple. We have two different met- 
als with junctions at temperatures T and T + e in two heat reservoirs. 
The point is that the temperature difference e 7^ generates by itself 
an electro-motive force. In other words, we can counterbalance by a 
non-zero electric field E and still have zero electric current. It can be 
measured in various ways, either by inserting condensor plates in one of 
the metals and seeing a potential difference or by inserting a battery in 
one of the metals and measure the voltage necessary to cancel the electric 
current. 

The simplest realization in our model is again the one of Section |21 
Instead of thinking of one wire (with M cells each consisting of two 
energy levels) we can think of two different wires (each with M cells 
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containing one energy level). What we have called the energy cur- 
rent in (j2.3p now becomes the electric or particle current. We have 
Pa = P,Pb = /3 + A. We insert an electric field of strength NE at the 
bond (iV, 1) (the -B— junction). The differential thermo-electric power 
or Seebeck-coefficient can now be computed by first solving for E the 
equation 

NE/3b + /3Be = /3e (6.15) 

which cancels the particle current and then putting v = —E/VT, cf. 
(j6.14|) . One easily finds that now 

v = -ep (6.16) 

6.2.2. Peltier effect. The Peltier effect is the production/absorption of 
heat at the junctions by the presence of an electric current across. Even 
without establishing a temperature difference, one generates a heat flow. 
The electric current is generated by an electric fleld or by an established 
electro-chemical gradient. To say it differently, to maintain a uniform and 
constant temperature throughout the system, absorption of heat from the 
reservoir is necessary. That is the Peltier effect used for example to cool 
an external device. 

Let us again look at the toy model of Section |21 By deflnition, tt/S is 
the isothermal entropy flux per unit electric current. Since the energy 
dissipation into the reservoir at inverse temperature /? equals e , we get 
71 P = e(3. Comparing with ()6.16|) . we see again ti (3 = —v, the second 
Thomson relation. 

In the general model particle current and heat current are not propor- 
tional and there is no simple analogue to (|6.15|) . Fortunately, the Onsager 
relation gives v = — vr [3 automatically. 

7. Beyond linear order, final remarks 

Despite their general appearance, fluctuation relations like ()5.9|) . ()5.11|) 
or ()6.7|) are not suflicient to build a systematic perturbation theory for 
the physical currents beyond linear order. The reason is that, for say 
the second order term in a current, we need to know the linear order 
for a time-symmetric observable (the square of the current) and that in- 
formation is absent from the fluctuation relations. There are however 
alternatives. Instead of going through general procedures, we outline 
here how for the model at hand one can compute higher order correc- 
tions to ()6.9p and ()6.12|) . The more systematic expansion will be part of 
another manuscript. 

As is clear form Section |3J all steady state averages of the form (J) (where 
J is some current), can be computed from the steady state average of the 
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particle current J^'^ over a bond (A;, /c + 1) in the bulk. That is directly 
clear from (j4.5p for the electric current J = but it also includes 
J = Jp the Peltier current as we can see from (j4.4|) . We thus only need 
to consider one bond {k, k + 1) and one energy level e & Sk and estimate 

We consider the case where we take a symmetric exclusion hopping pro- 
cess in the bulk of each line, i.e., Ek = for k ^ M,N. Due to the 
symmetric exclusion the occupation profile on energy level e is linear and 
the steady state particle current on that level is exactly given by the 
difference between the corresponding occupations at the two junctions in 
the same sense as . More precisely, for = 1 , . . . , M — 1 

, _ Prob[r/(e,l) = l]-Prob[r7(e,M) = l] 
M-1 

and similarly for the other bulk bonds. Hence all boils down to getting 
an expression for the occupations 

Prob[r/(e, A;) = 1], e e £k, k = I, M, M + 1, N 

To proceed we can take the same strategy as in ()2.Hj) - ()2.5|l and suppose 
that M, are sufficiently large so that the occupation numbers are given 
by the Fermi-Dirac equilibrium distributions. The driving parameters 
are present in these occupations through the different temperatures at 
the junctions and/or through the shifting of the energy levels by the 
presence of a local electric field. We do not continue giving the details 
of the computation but it should be clear that the computation is most 
pleasant in the grand-canonical ensemble where the occupations, and 
hence all currents can be straightforwardly computed to any given order. 
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